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Abstract 

The method of potential envelopes is used to analyse the bound-state spectrum 
of the Schrodinger Hamiltonian H = —A — v/{r + b) , where v and b are positive. 
We established simple formulas yielding upper and lower energy bounds for all the 
energy eigenvalues. 
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1. Introduction 



The cutoff Coulomb potential /(r) given by 



v/{r + b) 



(1.1) 



is an approximation to the potential due to a smeared charge distribution, rather 
than a point charge, and is appropriate for describing mesonic atoms [1]. Many 
authors have studied the eigenvalues Ene, n = 1, 2, 3, £= 0,1,2, ... generated by 
the cutoff Coulomb potential and have tried to estimate them. For example Ray 
and Mahata [2] applied the method of large- expansion to approximate the bound 
states energies from n = 1 to n = 4 . Mehta and Patil [1] rigorously analysed the 
S-wave bound-state eigenvalues of this potential as a function of b . 

In this paper we offer an elementary proof that the cutoff Coulomb potential 
has infinitely many discrete negative eigenvalues Ene,n = 1,2,3, = 0,1,2,... 
by using the comparison methods. We then use the comparison theorem and the 
envelope method [4-9] to obtain simple upper- and lower-bound formulas for all the 
eigenvalues. 

2. The discrete spectrum : Scaling 

The Hamiltonian for the problem is given by, 



A concern might be that, for suflBciently small coupling v , the potential, like a 
square well, might not have any discrete spectrum. However, the Coulomb tail averts 
this problem. It has been proved [5] by general methods that for any potential, 
like —v/{r + b), which is negative and decays at infinity slower than l/r^~^ , the 
corresponding Hamiltonian operator has infinitly many negative eigenvalues. The 
specific result for our problem may also be obtained by an elementary application 
of the comparison theorem, as we now show by the following argument. We note 
that the potential can be written 



H ^ -A- v/{r + b), V, b>0. 



(2.1) 



r^(r -|- b) 



(2.2) 



It therefore follows that 



V 



—V 



vb 
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and consequently 
where 

A = ((£+^)2 + ^6)^-l (2.4) 

Hence, we see that the effective potential associated with /(r) is bounded above 
and below by Hydrogenic effective potentials with discrete negative eigenvalues. 
This implies that the potential V has infinitely many negative discrete eigenvalues 
Ene bounded by 



These bounds are asymptotically close for large n. Another upper bound is provided 
by the linear potential since /(r) < — f + ^r. Hence, 

Eni<-l + {^)-^Sr,iil), (2.6) 

where Sniil) are the eigenvalues of the Hamiltonian —A + r for linear potential. 

For the S-states the radial equation may be transformed into Whittaker's equa- 
tion which has known exact solutions [3] . The general solution is written [3] in terms 
of the confluent hypergeometric functions M[x,y,z] and U[x,y,z] where, 

oo 

U{x,y,z)^^J e-'^H''-\l + t)y-''-^dt = z-''2Fo[x,l + x-y;;-l/z] (2.7) 



and M[x,y,z] = iFi[x]y; z]. Mehta and Patil [1] used the bounded property 
of the radial wave function and the boundary conditions to demonstrate that the 
eigenvalues are determined by the equation 

U[l-v/{2V^),2,2bV^] = 0. (2.8) 

As an alternative, we shall apply the envelope method to approximate all the eigen- 
values. We first reduce the complexity of the problem by the use of scaling argu- 
ments. If we denote the eigenvalues of H = —u)A — v/{r + b) by S{uj, v, b) , and 
consider a scale change of the form s — r/a , and choose the scale a = u>/v, then 
it is easy to show that, 

£{u;,v,b) = -S{l,l,-). (2.9) 
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Hence, the full problem is now reduced essentially to the simpler 1-parameter prob- 
lem 

H = -A-l/(r + b), £ = £{b),b>0. (2.10) 

3. Energy bounds by the Envelope Method 

The Comparison Theorem of quantum mechanics tells us that an ordering be- 
tween potentials implies a corresponding ordering of the eigenvalues. The 'envelope 
method' is based on this result and provides us with simple formulas for lower and 
upper bounds [6-9]. We need a solvable model which we can use as an envelope 
basis. The natural bases to use in the present context are the hydrogenic and linear 
potentials 

/j.(r) = sgn(q)r^, where q=— 1,1. (3.1) 

The spectrum generated by the potential h{r) may be represented exactly by the 
semi- classical expression 

Sniiv) = min{6' + vhniis)}, (3.2) 

s>0 

where the 'kinetic potential' hne{s) associated with the power-law potentials (3.1) 
are given by [9] 

h{s) = i2/q)\qSlf/{2 + g)|(''+^)/2s-'/^ (3.3) 

and is the exact eigenvalue of — A -|- sgn(q)r^, that is to say, corresponding 
to the pure-power potential with coupling 1. If we now consider a potential, such 
as /(r) , which is a smooth transformation /(r) = g{h{r)) of /i(r), then it follows 
that a useful approximation for the corresponding kinetic potential fniis) is given 

by 

fniis) ^ g{hnt{s)) (3.4) 

If g is convex in (3.4), we get [6-9] lower bounds ( ~=> ) for all n and £, and if g 
is concave we get upper bounds ( ~=< ) for all n and I . 

For the cutoff Coulomb potential, if we use the potential h = — 1/r as an 
envelope basis, then g is convex. An elementary calculation shows in this case that 



(3.5) 
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And if we use the potential h = r as an envelope basis, then g is concave, in fact 

—2v 

Thus in this application of the method we obtain upper energy bounds if we use 
h = — 1/r and lower energy bounds if we use h = r. The following remarks explain 
briefly how these results are obtained. 

We suppose for definiteness that the transformation g{h) is smooth and convex 
i.e g" > 0, then each tangent f^*\r) to g is an affine transformation of h 
satisfying 

/(*)(r) = a(t) + 6(t)%)</(r), (3.7) 
where the variables a{t) and b{t) are given by solving the contact equations 

f{t) = a{t)+b{t)h{t) (3.8) 

and f'it) = bit)h'{t), (3.9) 

which mean that the 'tangential potential', f^^\r) , and its derivative agree with 
/(r) at the point of contact, r — t. The potential inequality (3.7) implies the 
spectral inequality 

Ent[v) > va(t) + £ni{vb(t)). (3.10) 
The optimal lower bound thus obtained may then eventually [9] be re-written 

Er,i>miri{s + g{Ki{s))] . (3.11) 

In the complementary case where g is concave, the inequalities are reversed and one 
obtains upper bounds. 

For the power-law potentials h{r) = sgn(q)r'^ we can simplify (3.11) by chang- 
ing the minimization variable s to r defined in each case by the equation hniis) = 
h{r) so that g{h{r)) = f(r) = and the minimization (3.2), which yields eigen- 

value approximations for the Hamiltonian H = —u>A + f{r), where a; > 0, can be 
expressed in the form 

E^e ^ minL^ - -4-1 . (3.12) 



Spectral bounds for the cutoff Coulomb potential 



page 6 



We obtain a lower bound for Pni = Pni{—1) = (n + £), an upper bound for 
Pnc = Pntil), and a good approximation with the mean value Pni = P^ = 
|(-fn^(— 1) + -fn^(l))- These P-numbers are provided in Table (1). 

A natural question to ask is whether there exists a set of numbers {Pni} such 
that Eni = minj.>o 1^ + /(r)| exactly. We can see that the answer is "no" 
by an argument based on the 'concentration lemma' [10], which provides us with 
the relation between the concentration of the ground-state wave function and the 
size of the coupling constant v . More precisely, the wave function becomes more 
concentrated near the origin as v increases. Since for large values of the coupling v 
the "linear" upper bound (3.12) is very accurate (concentration near r = ), if there 
were one "exact" Pio , it would have to be the linear potential value Pio = Pio(l). 
But our upper bound is clearly above Ene for small values of v . Hence there are 
no such "exact" Pne- 

4. Results and conclusion 

We have derived a simple formula (3.12) for lower and upper bounds to the 
eigenvalues for the cutoff Coulomb potential. In Fig.(l) we plot the eigenvalue when 
(n, £) = (1,1) as a function of b for the case v = 1, accurate numerical values 
(dashed line), and our approximation with the average value Pn£ = |(Pn^(— 1) + 
Pne{l)) as stars. 

If we fix b and consider the Hamiltonian H = —A + vf{r), with eigenvalues 
S{v) , then from (3.12) we obtain the following explicit parametric equations for the 
corresponding approximate energy curve {v,£{v)}, namely 



2(Pn^) 

V — 



2 



■I ^ ' (4.1) 

These parametric equations yield upper bounds when Pne = Pnii^) lower bounds 
when Pne = {n + £), and a good approximation when we use the arithmetic average 
of Pne{—1) and Pne{^) ■ It is interesting, perhaps, that all these curves are scaled 
versions of any one of them; it is unknown if such a symmetry is true for the 
corresponding exact curves. In Fig. (2) we exhibit the graphs of the function E{v) for 
b = 1 along with accurate numerical data shown as a dashed curve. The main point 
of this work is to show that by elementary geometric reasoning one can obtain simple 
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semi-classical approximations for the eigenvalues. These results are complementary 
to purely numerical solutions and have the advantage that they are expressed simply 
and analytically and therefore allow one to explore the parameter space of the 
problem. 
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Table 1 The 'input' P-values P^^, P^^, and the mean = |(P„£(-l)+Pn^(l)) 
used in the general formula (3.12). 
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Figure 1. The eigenvalues £{b) of tlie Hamiltonian H = — |A — l/{r + b) for 
n = i = 1 (in atomic units h = m = 1). The continuous curves show the bounds 
given by formula (3.12), the dashed curve represents accurate numerical data, and 
the stars are the 'mean approximation' P^i = ^(Pni{—1) + Pnii^))- 
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Figure 2. Eigenvalue bounds (full-line) for the ground-state eigenvalue S{v) (n = 
l,i = )of the Hamiltonian H = —A + vf{r) (in units % = 2m = 1 ) for 6 = 1, 
together with accurate numerical data (dashed curve). The parametric equations 
(4.1) yield upper bounds when Pne = Pnei^) , lower bounds when Pne = Pn£{—1) 
and good approximation when Pni = i(P^£(— 1) + Pni{l)), shown as stars. 



